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Abstract 

The classical boundary-value problem of the Einstein field equations is studied with an arbi- 
trary cosmological constant, in the case of a compact (S 13 ) boundary given a biaxial Bianchi-IX 
positive-definite three-metric, specified by two radii (a, 6). For the simplest, four-ball, topology 
of the manifold with this boundary, the regular classical solutions are found within the fam- 
ily of Taub-NUT-(anti)de Sitter metrics with self-dual Weyl curvature. For arbitrary choice 
of positive radii (a,b), we find that there are three solutions for the infilling geometry of this 
type. We obtain exact solutions for them and for their Euclidean actions. The case of negative 
cosmological constant is investigated further. For reasonable squashing of the three-sphere, all 
three infilling solutions have real-valued actions which possess a "cusp catastrophe" structure 
with a non-self-intersecting "catastrophe manifold" implying that the dominant contribution 
comes from the unique real positive-definite solution on the ball. The positive-definite solution 
exists even for larger deformations of the three-sphere, as long as a certain inequality between a 
and b holds. The action of this solution is proportional to —a 3 for large a (~ b) and hence larger 
radii are favoured. The same boundary-value problem with more complicated interior topology 
containing a "bolt" is investigated in a forthcoming paper. 



1 Introduction 

In quantum gravity, as treated by the combined approaches of the Dirac canonical quantization 
and its dual, the Feynman path integral H, one studies (for example) the amplitude for n discon- 
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nected compact three-surfaces to have given spatial three- metrics, respectively h\p h 2 j, hfj. The 
amplitude is given formally by a path integral. In the simplest case n = 1, with only one connected 
three-surface, which could then be regarded as a spatial cross-section of a cosmological model, the 



amplitude is given by the 'no-boundary' or 'Hartle-Hawking' state [14| 

^Hu{hij) = j T>g^ u exp(-I E [g^]), (1.1) 

where the integral is over Riemannian four-geometries on compact manifolds-with-boundary, 
where the three-metric induced on the boundary agrees with the prescribed hij above. Here Ie 
denotes the Euclidean action [jl(| of the four-dimensional configuration, as in Eq. (|3.15| ) below. 

Naturally, in considering semi-classical approximations to this integral, one is led first to study 
the classical Riemannian boundary- value problem: namely, to find one (or more) Riemannian 
solutions to the classical Einstein field equations, possibly including a cosmological constant A, 
obeying 

Rfiv = Ag^, (1.2) 

which are smooth on the interior manifold, and which agree with the spatial three-metric hij at 
the boundary. Below, we shall see examples in which (different) such classical solutions exist, for 
a class of choices of the three- metric hij on a boundary diffeomorphic to the three-sphere 5 3 , on a 
manifold-with-boundary with the simplest possibility of a four-ball topology. 

Here, the boundary- value problem is studied within the class of biaxial Bianchi-IX models 
|7|, [DJ, which may be written locally in the form 

ds 2 = dr 2 + a 2 (r)( CTl 2 + a 2 2 ) + b 2 (r)a 3 2 . (1.3) 

Here <7j, i = 1,2, 3, denote the left-invariant one-forms on S 3 (see, for example, [fL3"l ). The boundary 
3- metric at a given value tq of r is then a Berger three-sphere |l9|, ^] , with intrinsic three- metric 

ds 2 = a 2 (r )(a 2 + a 2 ) + b 2 (r )a 2 . (1.4) 



Subject to condition fll.2| ), the metrics fll.3| ) are known in a closed form - they are the well-known 
family of Riemannian Taub-NUT-(anti)de Sitter metrics Q [j], [l3|]. For the topologically simplest 
case that we are studying in this paper, we insist that the S 3 with the given intrinsic metric Ql.4| ) 
bound a four-ball with smooth four-metric; this corresponds to requiring the Taub-NUT-(anti)de 
Sitter metrics to close with a regular "nut" and to have half-flat Weyl curvature. We find that the 
problem can be translated into an algebraic system of degree three which can be solved exactly 
to find the possible infilling geometries of this type and their actions for any boundary data a, b. 
Depending on a, b, one may therefore have three real roots, or one real root together with one 



complex conjugate pair, for this third-degree equation. A similar study was carried out in [17|, for 
large three-volume and small anisotropy, assuming a positive cosmological constant. However, we 
have been able to find all complex- and real-valued self-dual Taub-NUT-(anti)de Sitter solutions on 
the 4-ball and their actions for arbitrary values of (a, b) for both positive and negative cosmological 
constants. The case of zero cosmological constant can be obtained by taking the cosmological 
constant to zero (equivalently) from either the positive or the negative direction. 

Self-dual Riemannian Einstein spaces with a negative cosmological constant and of biaxial 



Bianchi-IX type have also been studied, for example, by Pedersen [19], in connection with the 
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conformal boundary- value problem with a three-sphere (topologically) at infinity. The more general 
case of generic (non-biaxial) Bianchi-IX models has been treated in this context by Hitchin [jn| 
and by Tod ||. 

Since a strictly negative cosmological constant A is needed to have any hope of incorporating 
gauge theories of matter with local supersymmetry in a four-dimensional field theory [||, ^4|], we 
have given a more detailed analysis for A < 0; the real and complex solutions have been classified 
completely in terms of the boundary values a, b, and the numerical behaviour of the solutions and 
Euclidean action I& has been investigated in greater detail. 



2 Einstein Biaxial Bianchi-IX Metrics of Riemannian Signature 



The general form of a biaxial Bianchi-IX four- metric is given by Eq. fll.3j ), Such metrics are invariant 
under the group action of SU(2) x U(l), whose Lie algebra is isomorphic to that of U(2). When 
one further imposes the Einstein equations with a A term, one arrives at the two-parameter Taub- 
NUT-(anti)de Sitter family § |, [l| : 

n 2 - T 2 dT 2 /\ 

ds 2 = P - dp 2 + , T „ (dip + cos 6 #) 2 + (p 2 - L 2 )(d6 2 + sin 2 6 # 2 ), (2.1) 
A ' p z — L z 



where 



A - ,r - 2Mp + L 2 + A ( L 4 + 2L 2 p 2 - ^p 4 ) . (2.2) 



Here L and M are the two parameters and < 6 < it, < < 2tt, < tp < Air/k (k is a 
natural number). When k = 1, the surfaces of constant p are topologically S 3 . The general form 
([2.1]), however, is only valid for a coordinate patch for which A / 0. In general A will have four 
roots. At the roots the metric degenerates to that of a round S 2 , and each such root therefore 
corresponds to a two-dimensional set of fixed points of the Killing vector field d/dtp. However, if a 
root occurs at p = \L\, the fixed-point set is zero-dimensional (as the two-sphere then collapses to 
a point). Such two- and zero-dimensional fixed point sets have been given the names "bolts" and 



"nuts" respectively [11]. The coordinate p ranges continuously from a root of A until it encounters 
another root of A, if there is any; otherwise p ranges from the root to infinity. In general, the bolts 
of the above Taub-NUT family of metrics are not regular points of the metric. For them to be 



regular, the metric has to close smoothly at the bolts, for which the condition is [18]: 

d ( A 



dpKp 2 -L 2 J^ pbM) 2kL- (2 - 3) 



For finite L / and k = 1, condition ( |2.3| ) leads to the self-dual Taub-NUT- (anti)de Sitter metric 
(see below) and the Taub-Bolt-(anti)de Sitter metric [Q, [l2|], which reduce to the Euclidean Taub- 
NUT |l5[ and the Taub-Bolt metrics []l8| respectively for A = 0. For the limiting cases L — > oo and 
L = 0, for which ( |2.1| ) is not well-defined, one can obtain regular solutions by suitable coordinate 
transformations and assigning correct periodicities to the coordinate parametrizing the S 1 fibre. 
For example, for A = 0, the Eguchi-Hanson metric [pfl and the Schwarzschild metric can be obtained 



from (|2,1|) in these limits [18|. In the next section we will encounter more examples of Bianchi-IX 



metrics that arise at these limits as we discuss the self-dual Taub-NUT-(anti)de Sitter solutions. 
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For a recent discussion on obtaining different biaxial Bianchi-IX metrics as various limits of ( |2.1| ) 
see §. 



2.1 Self-dual Weyl tensor and Regularity 

It will become evident below that the regular, positive-definite geometries of this type on a four-ball 
interior are Weyl half-flat, that is, having either self-dual or anti-self-dual Weyl curvature [^, ^] 
which requires: 

M = ±L(1 + -AL 2 ) (2.4) 
3 

For either sign, such metrics are often called (by physicists) half-flat - a name which we will be 
using in this paper. 

It is easy to check that, when the metric has a nut - or equivalently A has a (double) root at 
p = L (or p = —L) - then the two parameters L and M are related by 

M = ±L(1 + ^AL 2 ), (2.5) 

which is precisely the condition above for self-duality of the Weyl tensor. One can further check 
that Eq.fl2.3j) is automatically satisfied at the nut of the metric (2J). Therefore, within the Taub- 
NUT-(anti)de Sitter family, self-duality of the Weyl tensor will imply regularity at the origin. 

Without loss of generality, we will work with the positive sign of Eq.([0|) (and (2J3)). The 
condition ( |2,4| ) of half- flatness (or condition (2.5) of nut-regularity) implies that A simplifies: 

A = (p-L) 2 -±A(p + 3L)(p-L) 3 . (2.6) 
One can note that the terms involving A vanish as p — ► L, and that the metric near the origin is 



that of the A = case, namely Hawking's (Taub-NUT) solution [15 



„ 2 _ I P + L \ j 2 , AT 2 ( P~ L 



ds 2 = \J^Tl) dp + \p + LJ ^ + cos6d ^ + ~ L )( de + sin0 2 # 2 ) (2.7) 

Apart from the double root at p = L, A has two more roots or "bolts" at p = ±^(4L 2 + 3/A)- L. 
Beyond these two roots A is negative. So, the permissible range of p is from L to \J {AL 2 + 3 /A) — L 
(when L is positive) or from L to -y^L 2 + 3/A) - L (when L is negative). For a complete 
nonsingular metric, we have to check whether such points are regular, i.e., whether 



2{L- p)(p 2 + 3Lp + AL 2 )+3L\ _ 1 

(P=Pbolt) 



(P + L) 2 ) (o=0 ^ 2L (2i 



which is identically satisfied at the nut (at p = L) as we have seen already. However, for p 
yj (4L 2 + 3/ A) — L, which is positive definite and greater than L, this would require: 



-AL--V(4A2L2 + 3A) = — (2.9) 

which can only happen for L negative (= — |y6/A) and hence we find a contradiction. (In other 
words, this would give a regular bolt at p = |y6/A - on the "other side" of p = L which is 
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negative, and hence cannot be reached continuously starting from L.) The same argument applies 
for the other root at p = — \J (4L 2 + 3/A) — L. Therefore, for positive cosmological constant, the 
half-flat (equivalently, regular-nut) solution, has a singular bolt at finite p which exists and cannot 
be made regular for any finite value of the parameter L / 0. However, in the limit L — > 0, one 
can obtain the standard round metric on S by suitable coordinate transformations; this is regular 
everywhere and has a regular nut at each of the two poles. For the L — > oo limit, on the other 
hand, one can obtain the regular Fubini-Study metric on CP 2 . This is regular everywhere and has 
a regular nut at the origin and a regular bolt at infinity . 

The case with a negative cosmological constant is different. Writing A = —A (A is positive), 
one has 

A = (p - Lf + -\(p + 3L)(p - Lf. (2.10) 

A now has two roots at at p = L and two others at p = (4L 2 — 3A) — L, which are beyond 
the permissible range of p > L (L positive) and p < L (L negative). So, for a fixed negative A, 
the one-parameter family of half-flat solutions of the Taub-NUT-anti de Sitter type are necessarily 
regular for p G [L, oo) for any finite, non-zero value of L. Using similar coordinate transformations 
as in the case of positive cosmological constant, one can obtain the canonical metric on H 4 and 
the Bergman metric on CP 2 in the limits of L — > and L — > oo respectively, both of which are 
nut-regular at the origin. 

So, for both A > and A < 0, in this half-flat Taub-NUT-(anti) de Sitter family any hyper- 
surface of constant p (which is a Berger sphere with a given 3-metric) bounds a four-ball with a 
regular 4-metric. The singular bolt at p = iy 7 (4L 2 + 3/A) — L, in the case of A > 0, poses no 
problem as the Berger sphere lies between the (regular) nut and the (singular) bolt when we are 
interested in filling the Berger sphere with positive-definite regular metrics. It remains to see which 
classical 4-metrics of Taub-NUT-(anti)de Sitter type can be given on the four-ball inside a given 
Berger sphere. In other words, how many members in this one-parameter (L) family of metrics are 
there for which a given Berger sphere is a possible hyper surf ace at some constant pi 

In addressing this question we can include complex solutions in the interior of the 3-sphere. 
For this we let p and L and hence the 4-metric be in general complex-valued as long as the 4- 
metric induces the prescribed positive-definite 3-metric on the boundary and is non-singular in the 
interior.^ We will return to this issue in the next section. Real positive-definite infilling 4-metrics 
therefore form a subclass of the complex-valued solutions to this boundary-value problem and do 
not necessarily exist for arbitrary boundary data. However, it is possible to find the necessary and 
sufficient condition on the boundary data for the existence of real positive-definite solutions as we 
will see in Sections 4 and 5. 



3 Infilling Geometries and their Action 

Following the discussion above, the problem of finding classical infilling regular/ (anti) self-dual 
Einstein metrics of Taub-NUT-(anti)de Sitter type on the four-ball bounded by a typical Berger 
sphere with two radii (a, b) can be translated into solving the following two equations in p and L: 

a 2 -p 2 + L 2 = 0, (3.1) 
1 For a detailed discussion on such complex- valued metrics on real manifolds see EtJ and references therein. 
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and 

b 2 {p 2 _ L 2 ) - 4L 2 ((p - L) 2 - i A (p + 3L) (p - L) 3 ) = 0. (3.2) 

These are polynomial equations of degree two and six in the variables p and L and hence, by Bezout's 
Theorem (see, for example, would intersect at at most twelve points in (D 2 . However, trying 

to solve Eq.( |3.ip -(|3.2j) explicitly for p and L is not very easy and not to be recommended. However, 
there is a more "symmetric" approach, which will enable us to find explicit solutions, and first of 
all enables us to count the number of infilling geometries. 



3.1 Number of Infilling Geometries 

Theorem: For any boundary data (a, b), where a and b are the radii of two equal and one unequal 
axes of a Berger 3-sphere (squashed 3- sphere with two axes equal) respectively, there are precisely 
three (modulo "orientation") self-dual Taub-NUT-(anti)de Sitter geometries which contain the 3- 
sphere as their boundary!^ 

Proof: Note that the system of equations ( |3.1| )-( pr2| ) admits the discrete symmetry (p, L) <-> 
(— p, —L). (This is what we mean by "orientation"). Make the substitution: 

X = P + L (3.3) 

y = p-l. 



The problem now reduces to solving 
and 

u2_ y 



a 2 = xy (3.4) 



b 2 = JL( x -y) 2 (3-2Axy + Ay 2 ). (3.5) 

Note that this preserves the discrete symmetry, (p, L) <-> (— p, —L) in (x, y) <-> (—x, —y). Substitu- 
tion of x in Eq. Q3.5| ) now gives the univariate equation: 

Ay 6 + (3 - 4a 2 A)y 4 + (5a 4 A - 6a 2 )y 2 - (3a 2 6 2 - 3a 4 + 2 a 6 A) = (3.6) 

which is a cubic equation in y 2 . The six solutions for y therefore will appear in pairs of opposite 
signs. Since a 2 is positive, this implies that the six corresponding solutions of (x, y) are of the form 
(x 1 ,y 1 ), (x 2 ,y 2 ), (x 3 ,y 3 ) and (-x^-yi), (-x 2 ,-y2), (-£3,-2/3)- 

Since the set of solutions have the symmetry (xi,yi) — > (— Xi, —yi), by applying the transforma- 
tion (x,y) — ► (—x,—y), we would obtain no new solutions for (x,y) (hence for (p,L)) and would 
reproduce the same set. Therefore there are six points in (D 2 where the two polynomials meet, i.e., 
six solutions for (p,L) which are related by the reflection symmetry (p,L) «-> (— p, —L). Hence the 
number of geometries modulo orientation for any given boundary data (a, b) is three. □ 



2 As we are dealing with finite L / for which (2.1) is well-defined, self-dual metrics of biaxial Bianchi-IX type 
(admitting S^-foliations) arising as limiting cases (L — > and L — > 00) of the self-dual Taub-NUT-(anti)de Sitter - 
as discussed in Section 2 - are precluded from the considerations here. Depending on two radii of the Berger sphere, 
possible infilling solutions of such metrics can easily be included when one considers the larger class of (self-dual) 
Bianchi-IX type metrics on the four-ball. 
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3.2 Real/Complex Roots and Real Infilling Metrics 

It is more convenient to rewrite Eq.( |3.6| ) for z = y 2 in the form 

Az 3 + (3 - 4AA)z 2 + (5A 2 A - 6A)z - (SAB - 3A 2 + 2 A 3 A) = (3.7) 

where we have denoted (a 2 ,b 2 ) by (A,B). It is easy to see that the real solutions for (x,y) (and 
hence for (p, L)) come from the positive roots of Eq. (|3.7| ). Depending on the boundary data (A, B), 
the number of positive solutions for z can range from zero to three. However, a positive solution 
of z and, hence a real solution (p,L), would not necessarily correspond to a real positive-definite 
infilling geometry. For it to qualify as a regular real solution in the interior, the metric should 
be well-defined for p taking values on the real line from L to p. As we discussed in Section 2.1, 
this can happen if and only if p and L have the same sign and \p\ > \L\. The latter condition 
is automatically satisfied as a consequence of a 2 being positive irrespective of whether the signs 
of p and L are similar or not. If p and L have opposite signs, the metric will not be positive 
definite for values of p within the interval (L, —L) on the real line, and will be singular at p = —L. 
Therefore, it would not qualify as a real positive-definite solution. Such solutions for (p, L) should 
be considered to correspond to complex metrics where p is generally complex-valued in the interior 
and real-valued on the given S 3 boundary and at L. (One can choose contours for p in the complex 
plane which avoid the singularity on the real line at —L and ensure regularity at L; see discussions 
in (TtJ.) We will discuss the dependence of positivity (and complexity) of z on the boundary data 
(A, B) in detail for A < in Section 5 to explore the structure of (real) solutions and their actions 
and shall return to this issue. We will show that for a given boundary data (A, B) the real regular 
infilling solution is unique and exists when a certain inequality holds between A and B. 

For negative or complex roots of Eq.(|3.7|), the interior metrics are necessarily complex- valued. 

3.3 Explicit Solutions 

We now write down explicit solutions for our boundary- value problem. Rewrite Eq. (|3.7| ) with the 
conformal rescaling 

A\A\ -> A, 

B\A\ -> B, (3.8) 
z|A| -» z, 

so that A, B are now dimensionless. For positive cosmological constant, this gives 

f(z) := z 3 + (3- AA)z 2 - {5A 2 - 6A)z + (3AB - 3A 2 + 2A 3 ) = 0, (3.9) 

and, for negative cosmological constant, 

g(z) := z 3 - (44 + 3)z 2 + (5A 2 + 6A)z + (3AB - 3A 2 - 2A 3 ) = 0. (3.10) 

Since we will be more interested in the case of negative cosmological constant, we only write down 
the solutions of g(z) = 0. These are 

1 |4 2 + 4A + 1 4 
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and 



Z 2 = _JL5L_3 + 1 + _^ + _^/3 _3/Q_3 /g 12 x 

12 A yo 3 2 \6 V r 



z 3 = z%, (3.13) 

where 



Q = 72 A 2 + 8 A 3 + 216 A - 324 AB + 216 + 36 yj AB (81 AB - 108 - 108 A - 36 A 2 - 4 (3.14) 

The solutions of /(z) = can be found by just nipping the signs of A and B. It is important to 
note here that, the appearance of explicit factors of i in the expressions ( 3.11| )-( p.l3 ) for the roots 



may not be an accurate guide as to whether a root is real or complex. In Section 4, we will see how 
to analyse the roots of g(z) (in a way that would enable one to determine whether the roots are 
complex or real, for different values of (A, B)) without invoking the explicit solutions ( 3,11| )-( |3.13 ), 



However, we do need the explicit solutions in order to calculate the action for such geometries. 

3.4 Action Calculation 

The Euclidean action is given by |Uj] : 

8ttGI e = - I d A Xyfg i-R — A) - [ d 3 xVhK (3.15) 
JM V 2 / JdM 

where M is the Riemannian 4-manifold with boundary dM, and the 4-geometry agrees with 
the specified three-metric hij on the boundary, i.e., 9ij\ gM = hij. K is the trace of the extrinsic 
curvature Kij of the boundary. 

For a classical solution, obeying the Euler-Lagrange equations of the action ( 3.15| ) (the Einstein 
equations = Ag^ u ), one has 

8itGI e = -A d A x^fg~- [ d 3 xVh~K. (3.16) 
J M JdM 

Now J M d A Xy/g is just the volume of A4, which for the self-dual Taub-NUT-(anti)de Sitter metric 
gives 

- A / d^g = ~^ 2 L(p- L) 2 (p + 2L) = ir 2 [y 2 - A) [y 2 -3 A). (3.17) 
JM 6 6 

The trace of the extrinsic curvature is (2^- + \)/N for any metric of the form 

ds 2 = N(r) 2 dr 2 + o 2 (r)(<r 1 2 + a 2 2 ) + b 2 (r)a 3 2 (3.18) 

giving 

-/ d 3 x^hK = --(2- + b ^](l67T 2 a 2 b). (3.19) 
JdM N \ a bj 

For the self-dual Taub-NUT-(anti)de Sitter metric the surface term in Eq.( |3.15| ) is therefore 

8 o (y 2 -A) (Ay* + AAy 2 + 3y 2 + 9A-8AA 2 ) 
d 3 xVhK = -ir 2 ^ '-±-2- — '-. (3.20) 



/ 

JdM 



3 A 
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The total action is then: 

Using Eq.(p.7|), Ie can be further simplified to give: 

SttGIe = -vr 2 (7 A A 2 - 12 A + 12 z + 3Az 2 - lOAAz + 35) (3.22) 
3 

which in terms of the scaled variables gives, in the case A > 

Q 

87rG(AI E ) = -vr 2 (7A 2 -12yl + 12z + 3z 2 -10^z + 3 J B) (3.23) 
3 

In the case of a negative cosmological constant, this is just: 

8ttG(XIe) = --vr 2 (7^ 2 + 12,4 - 12z + 3z 2 - 10 Az -3B) (3.24) 
3 

Clearly, the substitutions (|3.3| ) have led to considerable simplifications. The action is a quadratic 
function of z - the solution of the third-degree equation - and hence would, in general, be complex 
(real) for z complex (real). Corresponding to every solution of Eq. (|3.9|) or Eq. fl3.10D , one can 
calculate the action merely by substitution in Eq.( 3.23| ) or Eq.( |3.24| ) - one does not need to find y 



(or x) at all to evaluate the action. As we discussed in Section 3.2, not all positive solutions for z 
correspond to real positive-definite solutions in the interior. Therefore all real- and complex-valued 
solutions in the interior corresponding to the positive solutions for z have real-valued Euclidean 
actions. This is because the two limits (here p and L) of the integral ( 3.17| ) of p lie on the real-line. 



Interestingly, for negative solutions for z, which correspond to p and L being purely imaginary, the 



actions are real as well. This is not surprising given that the middle expression of ( 3.17 ) is a quartic 



polynomial in p and L. Only the complex solutions for z would in general have complex actions. 

4 Vacuum Case 

Before embarking on solutions with non-zero cosmological constant, it is desirable to understand 
the A = case. This will also provide us with an example and illustrate the utility of the method 
developed so far. For this case, the regular Taub-NUT metrics are given by the Hawking solution 
(p. 7]). The corresponding (quadratic) polynomial equation and actions are obtained by simply 
setting A = in Eq.(|3,7|) and Eq.( |3.22j ) respectively. The two solutions and their actions are 

z = A±\[AB (4.1) 

and 

I E = ^(B±4VAB). (4.2) 



Following the discussions in Section 3.2, it is easy to check that the negative-sign solution of ( [4. 1| ) 
corresponds to a real positive-definite metric on the four-ball provided that 

a > b. (4.3) 
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This can be confirmed by solving for p and L directly (which is possible in this case). For the 
negative and positive signs of Eq. (|4.1[) they are 



2 v 'a? — ab' 2 

and 



,= ?i£4, L=U a ±^, (4.4) 



, L = - b J°±±R. (4 .5) 

2 Va 2 + a& 2 V a + b 

(Two other solutions are obtained by changing orientation.) Note that, in contrary to the statement 



made in [17|, positive-definite real infilling solutions do not exist for arbitrary boundary data 
although the actions ( |4.2| ) are real- valued for any value of A and B. The issue of positive-definite 
real solutions on the four-ball will be made clearer in the following sections where we discuss the 
case of non-zero cosmological constant (which includes A = as a special case). 



5 Solutions for A < and their Actions 

We now treat the case of a negative cosmological constant in greater detail and systematically 
discuss the structure of the solutions of 

g(z) := z 3 - (4A + 3)z 2 + {5A 2 + 6A)z + (3AB - 3A 2 - 2A 3 ) = (5.1) 

as functions of boundary data (A, B). Recall that (p,L) can have a real solution only when z is 
positive. Therefore, by studying the solutions of Eq.( ^.l|) , we are able to see the real and complex 
solutions of this boundary-value problem quite readily. However, before doing so we discuss the 
possibility of real positive-definite infilling solutions. 

5.1 Real, Positive-definite Infilling Solutions 

For a given boundary data (A,B), let the infilling solutions be denoted by (p,L). Recall from 
our discussion in Section 3.2 that a real positive-definite infilling solution exists if p and L have 
like signs (and if |p| > \L\ which holds automatically). Assuming without any loss of generality, 
that both p and L are positive, the equivalent requirement in the variables x, y is x > y. It then 
follows from Eq,(|3.4j) that the corresponding requirement on the possible solutions for z = y 2 is the 
following: 

z < a 2 . (5.2) 

Therefore, only positive roots of Eq. ( |3T9|) or Eq,( |3.10|) within the interval (0, A) give positive-definite 
real solutions on the four-ball which are regular everywhere including the origin, i.e., at p = L. 
Note that the requirement (|5.2| ) is independent of the size or sign of the cosmological constant. 

We now analyse the case of negative cosmological constant and find the necessary and sufficient 
conditions on the boundary data for roots of Eq. fl5.1| ) to satisfy the inequality (|5.2;). 
Lemma: There is a unique real, positive-definite self-dual Taub-NUT-anti de Sitter solution on 
the four-ball bounded by a given Berger-sphere of radii (a, b), if and only if 

b 2 < -a 2 (2a 2 \A\ + 3). (5.3) 
3 
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Proof: Recall that, for any function F(x), if F(a) and F(b) have unlike signs, then an odd number 
of roots of F(x) = lie between a and b (see, for example, ||). Depending on the boundary data 
there are two distinct possibilities for the sign of the constant term in Eq,(|5.l|). For B < ^^4(2^4+3), 
g(0) is negative and g(A) = SAB and is strictly positive. Therefore, for B < ^A(2A + 3), there 
will be either one root or three roots in the interval (0, A). By direct differentiation, one can show 
that g'(z) and g"(z) are strictly non-negative within the interval [0,A]. It therefore follows from 
Fourier's theorem that the number of roots in the interval (0, ^4) is one.f] If, on the other hand, 
B > ^A(2A + 3), there will be no root within the interval (0, ^4). □ 



Note that ( |5.3j ) gives (4.3) in the |A| = limit. There are certain particular features of the 
coefficients of g(z) that simplified our analysis. Firstly, it is only in the constant (z°) term that 
B appears. Also, and more importantly, except for the constant term the coefficients of g{z) are 
either positive or negative definite for arbitrary (A,B). (In the case of a positive A, we are not 
fortunate enough to have this particular advantage - the corresponding analysis is therefore just 
a step lengthier, although the methodology developed in this section equally applies. Possible 
positive-definite real solutions are unique as in the case of negative cosmological constant.) 



5.2 Structure of Roots 

We now discuss how the roots of Eq.(|5.1|) change as the boundary data is varied. Since for real 
(complex) roots the action is real (complex) this would enable us to study the real (complex) 
actions. In fact, as we will see below, regions of physical interest are covered by real solutions of 
Eq.Q. 

As is apparent from the previous discussions, we must look carefully at the two different cases, 
for which the constant term is positive or negative. The condition determining the reality of the 
roots of g(z) is following (see, for example, Q): 

^ < in" ^ vf^ (three real roots) 

B = M ^ (three real, at least two equal) (5-4) 

^ > w ( one rea ^ * wo com Pl ex roots) 

However, for the discussion below one further needs to know Descartes' Rule Q: the equation 
g(z) = can not have more positive roots than g(z) has changes of sign, or more negative roots 
than g{—z) has changes of sign. (It is also convenient to remember the elementary fact that the 
product of the roots of a polynomial equation of odd degree is minus the constant term.) 



Case 1: B < \A{?> + 2 A) 

In this case (3AB — 3A 2 — 2A 3 ) is negative and hence g(z) has three changes of sign and g(—z) 
has none. Therefore, there will be either (1) three positive roots of g(z) or (2) two complex and 
one positive roots, depending on which condition of ( |5.4[ ) (A, B) satisfies. One of the positive roots 

3 Fourier's theorem: If F(x) be a polynomial of degree n and fi, F%, F n are its successive derivatives, the 
number of real roots R which lie between two real numbers p and q (p < q) are such that R < N — N' , where N and 
N' (N > N') respectively denote the number of changes of sign in the sequence /i,/2, •■•/n, when x — p and when 
x — q. Also ((iV — N') — R) is an even number or zero. See, for example, H. 
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in the former case and the positive root in the latter case will correspond to real positive-definite 
regular solutions in the interior. 



Case 2: B > §A(3 + 2 A) 

In this case, (3AB — 3A 2 — 2A 3 ) is positive and hence g(z) has two changes of sign and g(—z) 
has one. So g(z) has either (1) two positive roots and one negative or (2) two complex and one 
negative, depending on which condition of (|5.4[) (A, B) satisfies. 



Case 3: B = §A(3 + 2 A) 

For this special case (3AB — 3 A 2 — 2 A 3 ) = and hence the valid roots are the solutions of 

z 2 - (4A + 3)z + (5A 2 + 6A) = (5.5) 

which will have either two positive or two complex roots. (The other root is z = which is not 
physically meaningful.) However, neither of the positive roots, which exist if A < |, will be less 
than A and hence there is no real infilling solution in the interior. 



We now combine Descartes' rule with the conditions of (|5.4|) (Fig. 1). We plot the curves 
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Figure 1: Structure of Roots (for A < 0) 



B = ^^4(3 + 2 A) (green line) and B = ^j- ^ A+ ^ (blue line). For physical interest, we have also 
plotted the line A = B which represents isotropy. The two curves divides the (B, A) plane in four 
regions. 
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Region I 

This is the region bounded between the two curves on the left-hand side of their intersection point, 
with the green line included. One root is negative and two others positive. The positive roots 
become equal on the blue line and the negative root become zero on the green line. None of the 
positive roots corresponds to a real infilling solution in the interior of the Berger sphere. 

Region II 

This is the region above the intersection point of the two curves. So there is one pair of complex 
roots and a negative root. On the green line there are only a pair of complex roots, the real one 
being zero. 

Region III 

Two roots are complex and one is real and positive. The positive root corresponds to a real 
positive-definite infilling solution in the interior. As above, on the green line, the real root is zero. 

Region IV 

This is the region bounded by the two lines below their point of intersection, with the blue line 
included. All roots are positive and one of them corresponds to a real infilling solution. Two roots 
are equal on the blue line. 

5.3 Numerical Study 

In the previous section we have studied the general behaviour of the roots systematically. Although 
B appears only in the constant term of g(z), its role is rather crucial in determining the sign and 
reality of the roots. In this section we will study the roots of g(z) and their corresponding actions 
numerically, and demonstrate an interesting connection with catastrophe theory. 

First we study the roots and actions for fixed values of A, as functions of B. From the analysis 
in the previous section (and from Fig. 1), one would expect to study three distinct generic cases, 
determined by the value of A less than, equal to and greater than its value where the two curves 
meet, i.e., at A = 3/2 (and B = 3). 

Case 1: A < 3/2 

For a fixed A < 3/2, by varying B one moves continuously from region IV (three positive roots) 
to the green line where one of the roots become zero (two still being positive), then into region I, 
where two roots continue to be positive, but the other is now negative, and then to the blue line 
where the positive roots become equal, and then to region II where they turn complex (the other 
still being negative). 

The behaviour of the three roots and the corresponding Euclidean actions are plotted in Figs. 
2 and 3 (as long as they remain real), as functions of B for A = 0.2 and A = 1 respectively. (The 
same colour is used to show the correspondence.) 
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Figure 2: Real Roots and corresponding Actions as functions of B {A = 0.2) 




Figure 3: Real Roots and corresponding Actions as functions of B [A = 1) 
Case 2: A = 3/2 

For the special value of A = 3/2, all roots are positive in region IV, where, on the point that the 
green and blue lines meet two become equal and the other becomes zero. On entering region II the 
two equal roots become a complex-conjugate pair and the other turns negative. The behaviour of 
the three roots is therefore similar to the previous behaviour, except that they all turn complex 
and become negative at the same value (B = 3) (Fig. 4). The corresponding actions also have a 
structure similar to those in A < 3/2. 




Figure 4: Real Roots and corresponding Actions as functions of B (A = 3/2) 
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Case 3: A > 3/2 

For A > 3/2 one starts in region IV with three positive roots, two of which become equal on the 
blue line and then turn complex and remain so in region II. The other root remains positive in 
region III and becomes zero on the green line to become negative in region II. See Fig. 5 (^4 = 3) 
and Fig. 6 (A = 100). 




Figure 5: Real Roots and corresponding Actions as functions of B when A > 3/2 (A = 3) 




Figure 6: Real Roots and corresponding actions as functions of B for large A(A = 100) 



The "Catastrophe Manifold" of I E 

One can see that certain patterns emerge for both g(z) and Ie- For fixed values of A, the real 
solutions of g(z) form a two-fold pattern as functions of B. One can check that this occurs for all A, 
large and small. The upper and lower folds turn over at B = and at B = respectively. 

This is easily understood with the help of Fig. 1: the curves (of two roots) in the upper fold meet 
when the two roots become equal at B = gj- ^ > i- e -' De f° re entering the combined region of 
II and III where they turn complex. On the other hand, on B = 0, g{z) has a double root equal 
to A (the other root is 2A) - therefore the lower folds turn over around the B = line. The 
surface g{z) = thus formed by placing such images successively (Fig. 11) is similar in structure 
to a cusp catastrophe manifold familiar in dynamical systems driven by a quartic potential with 
two control parameters (see, for example, |2(|). The minima of a quartic potential occur when its 
derivative (a polynomial of degree three) is set to zero and hence the catastrophe manifold represents 



15 



the equilibrium points of the system. The "catastrophe map" is then the part of projection of the 
catastrophe manifold onto the plane of the control variables bounded between the lines along which 
the folds turn over - in our case this is the region bounded between B = and B = ^j- ^ A+ j^ curves 
in Fig. 1, i.e., the union of region I and IV. However, one may wonder why a cusp does not appear 
in this case. This is because of the (particular) way in which A and B combine in the coefficients of 
g(z) and also because they are constrained to be positive - both facts are dictated by the physical 
configuration. One can obtain a cusp, however, by working with the new variable D(= AB) instead 
of B. A cusp will appear at A = —3, which is obviously outside the range of our physical interest. 

It is not difficult to see that Ie has the same catastrophe map, namely the union of region I and 
IV. However, the more important observation is that the relative orientation of the folds of Ie is 
the same as that of the corresponding roots of g(z), i.e., they do not cross each other. This pattern 
persists even when one gets very close to A = as shown in Fig. 7. For very small values of A the 
green and blue lines nearly overlap and coincide completely for A = 0. They therefore meet the red 
line at infinity. For higher values of A one gets a persistent behaviour as in Figs. 8-10. Therefore 
the "catastrophe manifold" of Ie is "diffeomorphic" to the one found above for z, i.e., the surface 
Ie is obtainable from that of z by a deformation which preserves the catastrophe map. This is not 
automatic or obvious given the form of Ie which is a quadratic function of z. Note, however, that 
the surface of Ie is not smooth where the upper fold occurs. 

The two-dimensional catastrophe map in a dynamical system with a quartic potential demar- 
cates the regions of three stable minima and one stable minimum. In our case they indicate a 
demarcation between the regions where there are three real Ie and one real Ie for the boundary- 
value problem. 




10 20 30 40 50 10 20 30 40 50 10 " 20 30 40 



Figure 7: Actions as functions of B (A = 25 x 1CT 6 , 25 x 1CT 4 , 3 x 1CT 2 ) 




Figure 8: Actions as functions of B (A = 3 x 1CT 1 , 8 x 1CT 1 , 1.2) 
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Figure 9: Actions as functions of B (A = 2.0, 5.0, 15.0) 




Figure 10: Actions as functions of B (A = 1500, 15000, 150000) 




Figure 11: The "catastrophe manifold" g(z) = : the surface of Ie is similar in shape and 
orientation except at the upper fold where it is not smooth. 

5.4 Large Radii and Small Anisotropy 

The observations that the space of Ie is diffeomorphic to that of z and that it does not intersect itself 
have immediate physical implications. This means that the dominant contribution will always come 
from one solution, namely the one represented by the green curve in Figs. 2-10 (and the meshed 
surface in Fig. 11). Note that this is the solution which gives a positive-definite infilling solution 
as long as b 2 < ^a 2 (2a 2 + 3). One can verify from Figs. 2-6 that this solution always takes values 
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within the interval (0,A). 

In cosmology, one is more interested in regions where b is not greatly different from a and 
usually when they are both large. All three roots are positive in this latter region. Taking A = B, 
the roots of g(z) are 

Zl = A + \-\ VT2AT9, 

z 2 = 2 A, ~ (5.6) 

z 3 = A + l + \ V12A + 9. 

It is easy to check that only z\ < A and, hence, z\ corresponds to the positive-definite real infilling 
solution in the interior. The dominant contribution to the path integral ( |1.1| ) will therefore come 
from z\ (which can also be checked explicitly in this case). The action of this solution is given by 

8irG(XI E ) = -^vr 2 (-9 + 3V12A + 9 + 4AV12A + 9) -^=vr 2 ^t, (5.7) 

3 V / y3 

so becoming more negative as A grows. Actions for Z2 and Z3 are positive and become more positive 
as A grows, as can be checked explicitly by direct substitution. 

6 Conclusion 

We have shown that, for a given boundary which is a Berger sphere (a squashed S 3 with two axes 
equal), there are in general three distinct ways in which one can fill in with a self-dual Taub-NUT- 
(anti)de Sitter metric. With suitable choice of variables, the problem of finding explicit solutions 
for such infilling metrics can be translated into a univariate algebraic equation of degree three and 
hence can be solved exactly. The Euclidean action Ie is a quadratic function of the solutions of 
this third-degree equation and hence, corresponding to every solution of this equation, the action 
can be found purely in terms of the boundary data - the two radii (a, b) of the Berger sphere, 
for both positive and negative cosmological constants. The positive and negative roots of this 
equation correspond to metrics for which actions are real - only complex solutions correspond to 
complex-valued actions in general. 

In the case of a negative cosmological constant, we have further discussed systematically the 
ranges of a and b where the solutions of the algebraic system lead to real- and complex-valued 
solutions in the interior and have studied the structure of the three roots. We found that one of the 
roots corresponds to a positive-definite infilling solution if and only if b 2 < |a 2 (2a 2 + 3). For small 
squashing, i.e., when a and b are of the same order and not too small, all three roots are positive. 
When a and b are exactly equal, this holds for small radius as well. We therefore investigated 
the roots and their corresponding actions numerically in this range (i.e., until two of them turn 
complex) as functions of the boundary variables a, b. We found that both the roots and the actions 
have structures similar to those of the cusp catastrophe in dynamical systems. The "catastrophe 
manifold" of Ie does not intersect itself which implies that the the dominant contribution will 
come from the positive-definite infilling solution. Further, the classical actions for large values of 
the radii (a, b) in their isotropic limit (of particular interest for cosmology) have been discussed. 
The (dominant) contribution coming from the positive-definite solution has an action proportional 
to —a 3 , thereby favouring large radii. 
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